RUB-TPII-03/06 



Fractional Analytic Perturbation Theory in Minkowski space and 
application to Higgs boson decay into a bb pair 

A. P. Bakulesfl and S. V. MikhailoJl] 

Bogoliubov Laboratory of Theoretical Physics, 
JINR, 141980 Dubna, Russia 



N. G. Stefanifl 

Institut fur Theoretische Physik II, Ruhr-Universitdt Bochum, D-4^180 Bochum, Germany 

(Dated: March 7, 2008) 

Abstract 

We work out and discuss the Minkowski version of Fractional Analytic Perturbation Theory 
(MFAPT) for QCD observables, recently developed and presented by us for the Euclidean region. 
The original analytic approach to QCD, initiated by Shirkov and Solovtsov, is summarized and 
relations to other proposals to achieve an analytic strong coupling are pointed out. The developed 
framework is applied to the Higgs boson decay into a bb pair, using recent results for the massless 
correlator of two quark scalar currents in the MS scheme. We present calculations for the decay 
width within MFAPT including those non-power-series contributions that correspond to the O (af)" 
terms, taking also into account evolution effects of the running coupling and the 6-quark-mass 
renormalization. Comparisons with previous results within standard QCD perturbation theory are 
performed and the differences are pointed out. The interplay between effects originating from the 
analyticity requirement and the analytic continuation from the spacelike to the timelike region and 
those due to the evolution of the heavy-quark mass is addressed, highlighting the differences from 
the conventional QCD perturbation theory. 

PACS numbers: ll.10.Hi, 11.15.Bt, 12.38.Bx, 12.38.Cy 



*Electronic address: bakulev@theor.j lnr.ru| 

^Electronic address: mikhs@theor.jinr.ru 

^Electronic address: stefanis@tp2.ruhr-uni-bochum.de 



1 



I. INTRODUCTION 



QCD perturbation theory in the spacelike (Euclidean) domain is based on a power-series 
expansion in terms of the running (effective) coupling a s (Q 2 ), (Q 2 = —q 2 > 0), which in 
one-loop order reads 

a,&) = ^ a(L) = - - L (1.1) 

2 

with b = 11 — Nf, L = ln(Q 2 /A 2 ), where A 2 = Aq CD , and with the "normalized" coupling 
3 

a(L) satisfying the renormalization-group equation 



da(L) 
dL 



-a 



(l + da + c 2 a 2 . . .) . (1.2) 



Here c\ = b\jb\ and c 2 = are auxiliary expansion parameters (see Appendix lAl and 

also Section III in Ref. [H). The one-loop solution of this equation suffers from an artificial 
singularity at L — 0, called the Landau pole. This prevents the application of perturbative 
QCD in the low-momentum spacelike regime with the effect that hadronic quantities, cal- 
culated at the partonic level in terms of a power-series expansion in the running coupling, 
are not well defined. Besides, from the theoretical point of view, such ghost singularities 
contradict causality rendering a spectral Kallen-Lehmann representation meaningless. On 
the other hand, in the timelike, Minkowski, region (q 2 > 0) the definition of the running 
coupling turns out to be difficult. The origin of the problem is that the QCD perturbative 
expansion cannot be defined in a direct way in this domain. 

Many efforts have been made since the early days of QCD to define an appropriate cou- 
pling parameter in Minkowski space in order to describe crucial timelike processes like the 
e + e~ annihilation into hadrons, quarkonium and r-lepton decays into hadrons, etc. Most 
of these attempts — see, for instance, 0, 0, 0] — were based on the analytic continuation of 
the strong coupling from the deep Euclidean region, where perturbative QCD calculations 
are safely performed, to the Minkowski space, where physical measurements are carried out. 
Over the years, it became clear that in the infrared (IR), the strong coupling may reach a 
stable fixed point and cease to increase. This behavior would imply that color forces may 
saturate at this low-momentum scale meaning that gluons decouple from quarks because 
they "see" them as a whole, i.e., in a quasi colorless configuration. Cornwal [EJ studied, 
within a vortex-condensate formalism, the formation of a mass gap, or effective gluon mass, 
that prevents the strong coupling becoming infinite at the Landau pole. Similar attempts 
were undertaken by other authors in subsequent years 0, 0, 0] using different techniques, but 
basing their arguments on the gluon acquiring an effective mass that works like an IR regu- 
lator in the low-momentum region. It was shown in j^] that this version of the IR-protected 
strong coupling can be related to the Sudakov factor for the non-emission of soft gluons, in 
the sense that gluons with wavelengths above some characteristic (nonperturbative) length 
scale, cannot resolve individual quarks because these segregate into a colorless mock-hadron 
state. 

In separate parallel developments, Radyushkin and Krasnikov and Pivovarov [ll[ 
have obtained analytic expressions for the one-loop running coupling (and its powers) directly 
in Minkowski space using an integral transformation from the spacelike to the timelike regime 
reverse to that for the Adler D-function (for more details, we refer the interested reader to 



12j, ll3|, ll4|). This sort of analytic coupling in the timelike region was rediscovered in the 



context of the resummation of fermion bubbles by Beneke and Braun [15| and also by Ball 
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the latter work in connection with techniques and applications 



Beneke, and Braun in 1 161 . 
to the t hadronic width. 

A systematic approach, termed Analytic Perturbation Theory (APT), has emerged in the 



last decade from studies initiated by Shirkov and Solovtsov [171 . Il8| . The main quantity of 



this framework is the spectral density with the aid of which an analytic running coupling is 
defined in the Euclidean region using a spectral representation. The same spectral density 
can be used to define the running coupling in the timelike region having recourse to the 
dispersion relation for the Adler function [ 



20 



These integral transformations, called R 
and D operations (see next section), coincide with those invented in [Io|, and provide 
the possibility to define simultaneously an analytic running coupling in both the Euclidean 
and the Minkowski space. Meanwhile this analytic approach has been extended beyond 

and important techniques for numerical calculations have been 
The approach has already been applied to the calculation 



the one-loot 



developed |22|, [23 



level J20J 

2. 
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25 



26, 27 



of several hadronic quantities, important examples being the inclusive decay of a r-lepton 
into hadrons [28|, 0, 3(|, the momentum-scale and scheme dependence of the Bjorken [3l[ 
and the Gross-Llewellyn Smith sum rule p2|, T decays into hadrons [33}, etc. Moreover, 
it has been extended to processes, like the 7*7 — > n transition form factor [i 



35 and the 



pion's electromagnetic form factor at the next-to-leading order (NLO) of QCD perturbation 
theory [ 



35, 36 



processes that contain more than a single perturbative scale, accounting 
also for Sudakov suppression. 

Overall, this analytic approach (see for some review [H, [37], H2]) does provide a quite 
reliable description of hadronic quantities in QCD, though there is also criticism [391 and 
alternative proposals and views of how to avert singularities in the running coupling |l3l. liol 

U H m, B, 0, 0, H S, 0, B B m, M — in particular concerning the deep IR region 

q2 <- whgj-g eventually the presence of a non- vanishing hadronic mass may become 



important [541 ] . It also suffers severe limitations in concept and application to processes 
beyond the leading order (LO) of QCD perturbation theory because it assumes that the 
only quantities that have to be analytic in the complex Q 2 plane are the running coupling 



56l . l57l . l58l . l59l . |60| that typical three-point 



and its integer powers. But it was shown in |5J 
functions, like the electromagnetic or pion-photon transition form factor, at the NLO level 
of perturbative QCD, and beyond, contain typical logarithms depending on an additional 
scale that serves as a factorization or evolution scale. These logarithms, though not affecting 
the Landau singularity, do contribute to the spectral density. This has led Karanikas and 
Stefanis (KS) [6l|, |62J to extend the concept of analyticity (the dispersion relations) from the 
level of the coupling and its powers to the level of QCD hadronic amplitudes as a whole. This 
generalized encompassing version of the analyticity requirement demands that all terms that 
may contribute to the spectral density, i.e., affect the discontinuity across the cut along the 
negative real axis — 00 < Q 2 < 0, must be included into the "analytization" procedure. The 
implementation of the KS analyticity postulate entails the extension of the original APT to 
non-integer (fractional) powers of the running coupling, which, in particular, encompasses 
the logarithms of the factorization (evolution) scale just mentioned. 

In this context it is worth emphasizing that fractional powers of the strong coupling were 
considered implicitly in J^. 1 The systematic development of the Fractional Analytic Per- 
turbation Theory (FAPT) for QCD in the Euclidean space was recently carried out in [l[ (see 
also [66] for a brief introduction) and was applied in [67J to the factorized part of the pion's 



1 It is interesting to recall here early attempts to study a spectral density amounting to fractional indices 
of the coupling in QED 64|. Such a spectral density was reinvented later within QCD by Oehme 65 1. 
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electromagnetic form factor that epitomizes three-point functions in perturbative QCD. The 
pivotal advantage of this scheme is the diminished sensitivity of the perturbative result on 
the factorization scale that parallels the strong renormalization-scheme independence, al- 
ready established within APT at the NLO level with respect to the same observable, in the 
exhaustive in-depth analysis of [36[ (an abridged version of which is given in j68| — see also 



[69l . I70l . I71j). The aim of the present investigation is to extend this analytic framework to 
the Minkowski space twining the two regions, spacelike and timelike, for any real index and 
any argument of the couplings and creating a new calculational paradigm for applications to 
hadronic observables in QCD perturbation theory. The formal discussion of the main charac- 
teristics of FAPT in Minkowski space (for which we use the abbreviation MFAPT) is carried 
out for two- and three-loop running-coupling parameters, investigating also the convergence 
properties of this type of expansion. To assess the consequences of MFAPT for observables 
and elaborate on its advantages in detail, it is best to study a quantity at a high-loop order 
of perturbation theory. To this end, we consider the correlator of two scalar bottom-quark 
currents, whose imaginary part Rs(s) is directly proportional to the decay width of a scalar 
Higgs boson to a bottom-antibottom pair. Considerable progress has been achieved with 
respect to this quantity during the last few years mainly thanks to the efforts of Chetyrkin 



and collaborators [72|, [73|, 174|. In the present investigation we will provide estimates for 
Rs{s) within multiloop MFAPT, using exclusively the MS scheme, and compare them with 
previous various results within conventional perturbative QCD up to the order O (af). The 
evolution effects due to the running of the strong coupling and the heavy-quark mass will 
be calculated up to three loops, borrowing the corresponding four-loop expansion coefficient 
for the Adler function from [74|. We emphasize in this context that our investigation is 
mainly meant to expose the conceptual advantages of the method, rather than to be used 
phenomenological tool for this quantity. 
The paper is organized as follows. In Sec. [Til we provide a mini review of the main features 
of APT, starting with the Euclidean region and completing the section with a discussion 
of the Minkowski domain. Though this material is mostly based on previously published 
works, its presentation here for any real coupling power is new and has not appeared before 
in the literature. Section II III is devoted to the generalization and extension of FAPT [l[ to 
timelike momenta, giving rise to MFAPT. In this section we present our main theoretical 
results and provide the reader with explicit (but approximate) two-loop expressions for the 
analytic powers of the timelike coupling. Similar formal expressions for any (higher) loop 
can be derived along these lines, with the three-loop case being outlined in Appendices |AJ 
IB"} and[C]). We close the discussion by giving analysis of the convergence properties of the 
perturbative expansion within our analytic scheme. Section IIVI contains an application of 
our framework on the correlator of two scalar currents of 6-quarks and the Higgs boson 
decay into a bottom-antibottom pair, presenting estimates for the width (actually for the 
quantity R$) of this process for different orders of the perturbative expansion, as specified 
above, and comparing our results with those obtained with the standard perturbative QCD 
expansion. Emphasis is put on the inherent advantage of our method to include into the 
analytic Minkowski couplings the crucial contributions stemming from the resummed 7r 2 
terms owing to the analytic continuation. Our conclusions are drawn in Sec. |V] where we 
also compile the benchmarks of FAPT in both the Euclidean and the Minkowski region. 
Important technical details are collected in four appendices. 
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II. CONCEPTUAL ESSENTIALS OF ANALYTIC PERTURBATION THEORY 



Here we introduce the main theoretical elements of APT, basing our considerations on 



371 . |75| |. the aim being to provide a comprehensive mini-review of the subject and equip 
the reader with all knowledge necessary for the extension to fractional powers in Minkowski 
space in Sec. IIHI As mentioned in the Introduction, the initial motivation to invent new 
couplings was the desire to interrelate the Adler D-function, 



calculable in the Euclidean domain, and the quantity R e 



a(e + e — > hadrons) 



a{e^e~ — > /i^/O 

R(s, /i 2 ) = £ r m {s/^) a"V) ^ R(s) = ]T r m a m {s) , (2.2) 

m m 

which is measured in the Minkowski region. Both quantities are considered in standard QCD 
perturbation theory, demanding that the couplings satisfy the renormalizat ion-group (RG) 
equation. In minimal-subtraction renormalization schemes the coefficients d n = d n (l), r m = 
r m (l), entering, respectively, the r.h.s. of Eqs. (12.1 p and (12.21) . are numerical constants. The 
functions D and R can be related to each other via a dispersion relation without any reference 
to perturbation theory. However, employing a perturbative expansion on the l.h.s. of Eqs. 
(12. ip and (12. 2p . one obtains, in fact, a relation between the powers of ln(s//Li 2 ) and ln(Q 2 /yU 2 ) 
in the coefficients r m (s//i 2 ) and d n (Q 2 /fi 2 ), while the powers of a s (/i 2 ) reveal themselves as 
numerical parameters. Upon setting fi 2 = Q 2 on the r.h.s. of Eq. (12.11) (or fi 2 = s in Eq. 
(12. 2p ). the coefficients d n (analogously r n ) become constants, whereas the coupling powers 
a™{Q 2 ) (equivalently, a™(s)) are part of the integral transformations (see below). But, if 
these coupling parameters are the standard running ones, then this connection fails at any 
loop order because of the Landau singularity in the Euclidean space. Questions arise whether 
analytic versions of both types of couplings may exist for which the above expressions could 
be connected. We shall see in the next step how this goal can, indeed, be achieved using 
non-power-series (functional) expansions within APT 12, 14, [HI, 0, 20, 37, [H, 75, 76l . 



The analytic images of the powers of the normalized running coupling, cf. Eq. (II. ip . in 
the Euclidean space can be defined by the formal linear operation Ae: 

Jo ° i* v 

where the spectral density is defined as 

f^((T) = -1m[a^(-(T)]. (2.4) 

The power (index) n denotes here only integer values, whereas the loop order is indicated by 
/ in parenthesis. We will show later that these relations are valid also for fractional powers 
(indices) v. 

Analogously, the analytic images of the normalized running coupling in Minkowski space 
are defined by means of another linear operation, Am, viz., 

A M [afj = 2l« with ^\s)= da . (2.5) 

J s ® 
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-s + ie 
• < 



%£ 



Q 2 -plane 



FIG. 1: Integration contour for the ii-operation in Eq. (|2.7p . 

These "analytization" operations can be represented by the following two integral transfor- 
mations: 

• D from the timelike region to the spacelike region 

-A® with 4P(Q 2 )=Q 2 r \^% da (2.6) 

Jo (a + Q 2 ) 



£>U„ 



and 



R for the inverse transformation (adopting the terminology of Shirkov, see, for instance, 



R[A®]=*$ with ^) {s )^±J~^^M d a, (2.7) 

where the last integral is evaluated along the contour shown in Fig. [TJ Note that these 
operations are connected to each other by the relation 

DR = RD = 1 , (2.8) 

valid for the whole set {A n , 2l n } and at any loop order. 

The operations A E and A M , which define, respectively, the analytic running couplings in 
the Euclidean (spacelike) and in the Minkowski (timelike) region are displayed graphically in 
Fig. [2j The logic of "analytization" enables a similar outcome with respect to the expansion 
of QCD amplitudes (depending on a single momentum scale Q 2 ) and their continuation from 
the Euclidean to the Minkowski space. As an example, consider the Adler function D on 
the r.h.s. of Eq. (12. II) . which is expanded in terms of a™(Q 2 ). The operation Ae, applied to 
D(Q 2 ) along the right arrow in Fig. [2]^a), maps it on a non-power-series expansion [l2|, 75| 
in the Euclidean region, termed Va, i.e., 

D(Q 2 ) = d n a n {l) {Q 2 ) A E [D] = V A with V A (Q 2 ) = £ d n A^(Q 2 ) . (2.9) 

n n 

Subsequently, one can apply the R operation, given by Eq. (12. 7p (bottom line in Fig. [2](a)), 
to obtain the quantity R in the Minkowski region: 

R [D A ] = K with K{s) = ^ d n Ql n (s) . (2.10) 
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FIG. 2: Implementation of analyticity in APT (a) and in FAPT (b) and (c). The index n in 
APT is restricted to integer values only, while in FAPT v can assume any real value, enabling the 
"analytization" of expressions like those shown in (c) and presented in Appendix [Cj 



On the other hand, the same expression for 7Z(s) (in Minkowski space) can be obtained 
following the left arrow in Fig. 12(a) by making use of the Am operation. This leads to the 
same analytic image of R(s): 



D(Q 2 ) = J2 d n a n {l) {Q 2 ) A M [D] = TZ with K{s) = ^ d r , 



Onto 



(2.11) 



The generalized concept of imposing analyticity to the QCD amplitude as a whole [61|, 162 



allows us to perform the "analytization" of any fractional (real) power of the coupling and — 
even more important — invoke the "analytization" concept on more complicated expressions 
that contain products of powers of the strong coupling times logarithms of a second pertur- 
bative scale, like the factorization or evolution scale — see for an illustration part (c) of Fig. 
[2] and for details Appendix O Then, all of the previous results, exposed via Eqs. (I2.3p -( l2~8lb 
can be generalized to hold for any fractional index (power) u, giving rise to the vector spaces 
{A u }, {SlJ that possess the property of index differentiation. 

Let us now turn our attention to the spectral density. At the one-loop level, we can derive 
by a straightforward calculation, based on Eqs. (11. ip and (12. 4p . a closed-form expression for 
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the spectral density; namely, 1 

= ~ ; > <P. = ™ I ) , L CT = In (a/A 2 ) . (2.12) 




Then, the one-loop couplings A\ L> and can be derived by substituting p± into Eqs. (12.31) 
and (ESD to get 

4 1} (Q 2 ) = I-^3T> ( 2 - 13 ) 



and 



a^fs) = - arccos [ ; Ls ) (2.14) 



with 

L = ln(Q 2 /A 2 ) , L s = ln(s/A 2 ) . (2.15) 

From these equations we infer that "analytization" (Ae) in the Euclidean case amounts 
to the subtraction (at the one-loop level) of the Landau pole, whereas in the Minkowski 
space the analogous operation (Am) means summation of 7r 2 -terms in all orders of the ex- 
pansion. To generate two-loop expressions for the analytic couplings, one can make use of 
the Lambert function, as shown by Magradze in [22|]. Still higher loops can be obtained via 
an approximate form of the spectral density and numerical integration (75| . It follows from 
this summarized exposition that the difficulties encountered with ghost singularities in the 
Euclidean strong coupling can be eliminated on account of causality ( "spectrality" [75j) and 
RG invariance. Hence, from the point of view of the analytic approach, the Landau pole re- 
mover (and analogously the compensation of singularities in higher loops) is not introduced 
by hand but ensues naturally as a corollary within the formalism without appealing to any 
nonperturbative physics as the origin of power corrections. Nevertheless, one may include 
into the spectral density power corrections of the form (M 2 /Q 2 )™, with M being, for exam- 
ple, a constituent quark mass, in an attempt to incorporate this way some nonperturbative 
effects. Such "analytization" approaches, following, however, different incentives, have been 
proposed in 5l|, |52J . 



III. MINKOWSKI VERSION OF FRACTIONAL ANALYTIC PERTURBATION 
THEORY 



Before studying the detailed procedure of the analytic continuation to Minkowski space, 
let us first make some important remarks about the strategy on how to generalize the ap- 
proach in order to include non-integer (fractional) indices v. Appealing to our detailed 
discussion in [l|, we note that it is not obvious how the direct way for the analytic continu- 
ation of the coupling powers in terms of Eq. (12. 3p can provide explicit analytic expressions. 
Therefore, to achieve this goal we have used instead in [l] another method, based on the 



2 An analogous expression for QED can be found in [6. 
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Laplace representation, which will be exposed in the next subsection. However, the situation 
in Minkowski space is different because of the absence of ghost singularities. In that case 
it turns out to be possible to employ the dispersion-relation techniques (cf. Eq. (12.51) ) in 
order to obtain analytic expressions in the timelike regime that are valid for any real index 
v. Indeed, using Eq. (I2.12p . we find in one- loop order 



"°° Art 1 r°° sin 



'arccos \ L/ a/ [L 2 + tt 2 )^ 



o ' v ^' * Jl~ (7f 2 + L 2 ) u/2 

This integral can be evaluated explicitly and provides the result 



(3.r 



sin 

rt!Hs) = — 



— 1) arccos \L S / a/ (L 2 + 



7r \u — 1) (L z s + n 2 y " 

that is completely determined by elementary functions [j], S]- Taking the limit v — > 1 in 
the above equation, one readily obtains Eq. (12.141) for 21] , while taking v = one finds 
2t« = 1. 

Let us consider now the spectral density pP{o~) beyond the leading-order approximation. 
At the /-loop level, pi (c) can always be presented in the same form as for the leading-order 
one, given by Eq. (12. 12f) . i.e., 

-Im [a {l) (-a)\ = v 

keeping, however, in mind that the phase and the radial part Ry) acquire now a multi- 
loop content. Suffice it to mention here that an explicit two- loop expression for the spectral 
density is derived in Appendix [Bj notably, Eq. ( IB 71) . In the same appendix we show that this 
expression, though approximate, is very close to the exact, but numerical, one. To be more 
specific, one should, strictly speaking, deal with the imaginary part of the appropriate branch 
of the Lambert function W_\ (see [22|) owing to the fact that the exact solution of the two- 
loop RG equation (given in Appendix |A] by Eqs. (1A4|) and (1A6j) ) can be expressed in terms 
of this function. To complete the exposed procedure, one should substitute the displayed 
spectral density into Eq. (12.51) and perform the integration. Recently, Magradze jH, 2?J has 



u)f n 7 \ut M sm l u ni)K a to ox 



(2) i2) 

published closed-form expressions for 21]/, 21;/ at the two- loop level by means of the W\ 
Lambert function. The dark side of this latter procedure is that it does not lend itself to an 
analytic evaluation of explicit expressions for 2lJ> 2 ' ) (L s ) for fractional indices, but yields (after 
integration) only numerical values. Beyond the two-loop level, explicit results are difficult 
to obtain. Numerical values of the quantities A n and 2l n for n = 1, 2, 3 at the three- loop 



level were given by Kourashev and Magradze in [24|. Very recently, Shirkov and Zayakin 
[33| have constructed a simple one-parameter model to emulate the first three (n = 1,2,3) 
analytic couplings at the three-loop level, both in the Euclidean and the Minkowski region, 
that claims an acceptable accuracy for practical purposes. 



A. Simultaneous derivation of FAPT in the Euclidean and Minkowski regions at 
the one-loop level 

In this subsection, we consider timelike and spacelike couplings in mutual comparison, ex- 
clusively in the one- loop approximation, omitting for this reason the loop label (/). The gen- 
eralization to higher loops will be presented in Subsection IIII CI Let us start the derivation 
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of the (M)FAPT analytic couplings {A u (k), % v (k)} with k being a logarithm of a momentum 
(or energy) scale by employing the Laplace- representation approach of It is useful to 
recall at this point that the initial sets {Ail, {2l n } have been constructed for integer values 



of the index and constitute vector spaces []J. Being able to create the elements A u , 2l„ for 
any real z/, one can complete the vector spaces {A}, making it possible to apply other 

linear operations to these spaces, e.g., differentiation with respect to the index v. 

It is important to appreciate that the generalization of all APT couplings to fractional 
(real) values can be performed within a single mould. To this end, we apply a differential- 
type RG-equation, like (11.21) . and, following [H, [75|, we first write 




(n-1)! 



d_ 

dk 



n-l 




(3.4) 



where the evaluation of the couplings a n (standard, n: power) and A n (analytic, n: index) 
in the spacelike region proceeds with k = L = ln(Q 2 /A 2 ), while their counterparts in the 
timelike region, 2l n , are calculated with the aid of k = L s = ln(s/A 2 ). 

To facilitate the transition to fractional index values, it is instrumental to employ the 
Laplace representation of both types of couplings — the analytic, A n {l), 2l n (Z), and the con- 
ventional ones, a n (L), — and define for k > 



a n (k) 
A n (k) 

2ln(fc) 





(3.5) 



To derive the last equation, we have used in Eq. (13.41) the one-loop RG equation. The key 
element in converting the APT couplings to the set {A n {k), 2l n (fc)}, valid for any index u, 
is the relation 

a»(t) \ j / ai(t) 

Mt) d =^ Mt) 



(3.6) 



that generalizes Eq. (13. 5p . From this, it is evident that a\{t) = 1. The explicit expression 
for A±(t) — worked out before in |l| — and that for the timelike coupling 2li(t) can be written 
as follows 

1 



( ax(t) 



At(t) 



t 



v-l 



r M 



\ sin(7rt)/(7rt) 



(3.7) 



/ 



Note that the first term in the Euclidean analytic coupling (second line in the above equa- 
tions) stems from the usual QCD term 1/L, whereas the 5-function term is related to the 
Landau-pole remover (second term in Eq. (I2.13P ). To reveal the particular features of the 
Laplace-conjugate images of these couplings, we show them graphically in Fig. [31 One sees 
from this figure that the Laplace conjugate of the conventional (normalized) coupling corre- 
sponds to a straight line at unity, while the analogous expression for the Euclidean coupling 
is represented by a Dirac comb (blue line) and the Minkowski one is a smooth and oscillating 
function (red line) dying out with t. 

Though we have initially assumed that k > and v > 0, these Laplace conjugates 
generate, in turn, the following images that can be analytically extended to any real u, L, 
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FIG. 3: Illustration of the Laplace images of the one-loop analytic couplings in the Euclidean 
(Ai(t)) and in the Minkowski space (2li(t)). 

and L a (we display again the loop label (1 = 1) explicitly for the sake of comparison later 
on) 

a (D = 7V > (3-8) 
\v — 1) arccos (l s / a/ (Xf + vr" 

7T (l/ - 1) (L2 + 7T 2 

where the Dirac comb gives rise to the transcendental Lerch function, F(z,u) (78[, that 
serves as a Landau-pole remover for any u, and the oscillating curve amounts to elementary 
functions, confirming the result given in Eq. (I3.2p . The last expression is a new result of the 
present analysis and has been derived by taking recourse to the last entry of Eq. (13. 7p . It 
is worth emphasizing that Av(L) and 2l£, (L s ) are entire functions of their corresponding 
arguments. 



sm 



al 1, (^.) = — ^ 7777± 2 ^- 1)/2 '\ (3-10) 



B. Properties of timelike vs. spacelike couplings for any real index v 



1. From inspection of the relations (12. 6p . (12.71) . and recalling that asymptotically as 
L — > oo both analytic couplings, A V (L) and 2lj, (L s ), tend to the same standard coupling 
a u (L), one may ask about the mutual behavior of these couplings for finite values of their 
arguments. Despite the asymptotic symmetry of the couplings, for finite arguments this 
symmetry is distorted 19, 2(j, albeit the couplings are equal at the origin, i.e., for Q 2 = s —>■ 
0, or equivalently, L = L s — » — oo. This "distorting-mirror" effect mm is an interesting 
property of the analytic couplings and was originally established for integer powers of the 
coupling. It is symbolically expressed through the operation 



(3.11) 
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(see Eq. (12.61) ). which we now generalize to be valid for any real index v. Its content may 
become evident from the following expressions for L = L s = 0: 

^(0) = Sm / (Z/ ~! )7r / 2] : Sti(0) = ^ 912(0) = 4 a 3 (0) = 0, 2U(0) = --^, (3.12) 

M0) = ~^W"- Mo) = l*Mo) = ^ Mo) = o, Mo) = -^q (3-i4) 

^(°) = ^^(°) = S^i (0) = -il' (3 - 15) 

where ((u) is the Riemann ( function. These couplings are interrelated by the equation 

v-i)C(y) 



A(o) 



tu~l 



21,(0) (3.16) 



with the coefficient in the bracket providing a quantitative measure for the magnitude of 
the distortion for any i/£l. For a graphic illustration of the "distorted mirror symmetry" 
effect, we refer the interested reader to [3, 19, 



2. As we have shown in [l[ for the Euclidean couplings, the parameters A- v play the role 
of the "inverse powers" of A\ that may be considered as the images of a~ v . This property 
extends also to the Minkowski region, so that the set {A- v , 21-,} corresponds to analytic 
images of a~ v in the Euclidean and the Minkowski space, respectively, for arbitrary v values. 
This allows us to demonstrate the "distorted mirror" effect in analytic form as follows 

2l (L s ) = l, 2U(L S ) = L S , 21„ 2 (L s ) = L2-tt 2 /3, 2L 3 (£.) = L S (L 2 S - tt 2 ) , . . . , , 
ML) = l, A-t(L)=L, A- 2 (L) = L\ A- 3 (L)=L 3 ,... [6 ' U) 

Note that the expressions for 2l_ n , given by Eq. (13.171) . can be linked to A~ n , in analogy to 
(13. lip , by means of the transformation 

A_ n (Q 2 ) = D[^M ■ ( 3 - 18 ) 

[See Eq. (19) in [73| for an earlier implicit derivation of these expressions, employing trans- 
formation (12.61) in terms of the powers of L s (or in terms of 2l_ n (s) in our notation)]. A 
useful all-order formula to analytically continue logarithms under the D transformation (see 
item 3 in Appendix [B]) has been presented in [63j . 

3. Moreover, the analytic couplings in both regions (Q 2 , s), or equivalently, (L, L s ), have 
the following symmetry (respectively, asymptotic) properties: 

2l m (L s ) = (-1)-2U-L S ) , An(L) = (-l) m A m (-L) for m > 2 , m e N ; (3 1Q) 
2l m (-oo) = Am(-oo) = 6 m>1 , 2l m (oo) = A m (oo) = for m E N . 

To reveal the details of behavior of the generalized analytic couplings 2l„(L s ) and A V {L) 
and make the above statements more transparent, we illustrate them in Fig. H] in terms of 
two graphics, which display the rate of change of these functions with respect to the index 
v and the argument L. Inspection of this figure in conjunction with Eq. (13.191) provides also 
information about how the zeros of the couplings occur. 
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-15 -10 -5 5 10 15 -15 -10 -5 5 10 15 



FIG. 4: Comparison of the Minkowski (left panel) and the Euclidean (right panel) analytic cou- 
plings, %l u (k = L s ) and A u (k = L), respectively, for incremental changes of the index v in the range 
2 to 3. 



C. Extension to higher loops and convergence of FAPT in the Minkowski vs. the 
Euclidean region 

The last topic of this section is to extend our results to higher loops and to discuss the 
convergence properties of FAPT in the timelike region. In the following exposition, we shall 
employ, for the sake of simplicity, a special notation for the derivatives with respect to the 
index v e M of the non-power expansion and define 

In our previous paper [l[ we have obtained an expansion of the two-loop analytic coupling 

Au\l) in terms of the one-loop analytic coupling A^P(L) (see Eq. (3.29) in [1]). Due to the 

linearity of this expression in we can immediately rewrite it for timelike couplings — 

see Appendix EJ Eqs. (lC2ap . (lC2bl) — by virtue of Eq. (jCip to obtain the two- loop result 

f(2) _ oiM i r, ,,-nofC 1 ) , „Hn2 , n_il o^W 



*W = aW + c x vV^l 2 + c\ [V 2 + V - 1] a£ 

2tE 4 + O (V^l 5 ) , (3.21) 



+ c? 



5 1 

V 3 + -V 2 -IV 

2 2 



where we employed the auxiliary expansion parameter c\ = 6i/6q. 

Next, we test the quality of the two-loop expansion of the analytic-coupling images of ap) 

and (ct(2)) 2 in the Minkowski region in comparison with the Euclidean one (refraining from 
displaying the latter because it is completely analogous — see Appendix [Cj. In doing so, we 
define the following quantities for v — 1 and v — 2: 



NNLO, i.e., retaining terms up to order cf 
N 3 LO, i.e., retaining terms up to order cf 



4 (v 3 + lv 2 -IV -|)st (1) 



A 4 (Sir) = A 3 + 2 m ; (3.23) 
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N 4 LO, i.e., retaining terms up to order cf (cf. Eq. (I3.26c[) ) 

2^ + 2 ci V St^g + c\ (3V 2 + 2V-2) 21^ 



^(2);MFAPT 



+ c\ (AV 3 + 7V 2 -6D - 1) 2£ 



(1) 



+ C? ( 5P 4 + y V 3 



8V 2 - 11V 



10 

y 



21 



(1) 

^=6 



(3.24) 



with analogous expressions for A 3 (A u ) and A 4 (A u ), obtained by using the evident substi- 
tution 21 -» A 
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FIG. 5: (a): The dashed red line corresponds to A 3 (_4i), whereas the solid blue line represents A4(„4i). 
(b): The dashed red line denotes As(2li) and the solid blue line A^Sli). The dotted black line indicates 
A [t ~ 2 (Ai) (left panel) and A 1 *" 2 ^) (right panel). 



Figure [5]illustrates the convergence quality of the FAPT expansion in the spacelike and the 
timelike regions in terms of the quantities A n , defined above, that encapsulate the deviation 
of our approximations from the exact results. In addition, we consider the auxiliary quantity 
A lt_2 (dotted black line) for the spacelike (L) and the timelike (L s ) regions, defined by the 
respective expressions 



(2)it-2 



(2) 



A (2lj,) 



2l< 2) 



(3.25) 



which, as one sees, result from replacing the exact spectral density (cf. (1B8I) ) by its second 
iteration (details are relegated to Appendix IB"1) . 

From Fig. [5] one observes that the convergence of Af ] and 2tS 2) at the N 3 LO of the 
expansion in the auxiliary parameter c\ is sufficiently accurate. In the Euclidean case (left 
panel), the errors — defined by A FAPT (*4i), (dashed red line) and A FAPT (Ai) (solid blue 
line) — induced by the non-power-series expansion are by a factor of two less than those in 
the Minkowski case (right panel): A^ FAPT (2li) (dashed red line) and Af FAPT (2li) (solid 



blue line). The largest error results around L 



0, but already for \L\ > 2 the 



uncertainty is less than a few per mil, rendering the convergence of the expansion highly 
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reliable. Deriving the two- loop explicit expression 



[v+l)(u + 2 



7/ 



Z/+1 

2! 

2z/ + 3 



Z> + Z> — 1 



21 



(i) 

u+2 



V A + V z + 1 + i/)£> 

3! 2 V ; 2 



21 



(i) 

i/+3 



(3.26a) 
(3.26b) 



^+1)^ + 2)^ + 3)^ , 3z/ 2 + 12i/+ll p3 _ z/ 2 + 2z/^ 2 



4! 

3z/ + 5 3z/ + 4 

V H 

2 6 



P 4 + 



P^ 



(3.26c) 



that is extended to the three-loop order of the running coupling in Appendix O (Eqs. (jC2al) 
to (lC2cj) ). a similar quality of convergence can be established for any desired fractional value 
of the index v. 
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FIG. 6: (a): The dashed red line corresponds to A^' FAPT (L), computed analytically via Eqs. (|C2a[) - (|C2cp 
(for C2 = 0), whereas the solid green line represents the exact expression for A^(L), cf. Eqs. (|2.3p and ()B8|) . 
(b): The dashed red line denotes Stj ' MFAPT (^s)j computed analytically via Eq. (|3.24l) . whereas the solid 
green line represents the exact result for %^p(L s ), given by Eq. (|B11[) . On both panels, we also indicate by 



* the results of Magradze, presented in 25l. |26||. 



In Figure [6] we show the results of the exact and the approximate calculation of 
A^' FAFT (L) and 2l2^' MFAPT (L s ). This figure also includes the results (denoted by the symbol 
*) obtained before by Magradze j25|. One observes that the agreement between his numer- 
ical estimates and our more elaborated calculations is excellent. We take the opportunity 
to remark that our opposite statements in notably Figure 4, were incorrect owing to an 
error in our code. This bug has now been eliminated, so that Magradze's numerical results 

(2) 

are fully supported by our calculation of the exact expressions for both analytic images A 2 
and 2lf . 

In concluding this section, it is worth listing the advantages of our calculation: (i) a 
good convergence of the non-power-series expansion, (ii) full control over the region in L s 
(correspondingly L), in which the maximum uncertainty occurs, making possible system- 
atic improvements, and (iii) and most important for practical applications, a considerably 
reduced uncertainty level of the expansion in the physically interesting region, say, beyond 
1 GeV, (which corresponds to L = 2) of only a few per mil. 



15 



IV. SCALAR CORRELATOR AND HIGGS BOSON DECAY INTO HADRONS 
IN MFAPT 



To bring out the contrast with the original APT and effect the advantages and the utility 
of the (M)FAPT machinery with regard to the conventional perturbative expansion, we 
consider in this section the decay width of the Higgs boson into a bottom-antibottom pair, 
T(H — ► bb), using a multi-loop approximation. We stress in this context again that, though 
there are no ghost singularities in the Minkowski space owing to the strong coupling, the 
analytic continuation from the spacelike to the timelike region entails so-called 'kinematical' 
7r 2 terms (see, for instance, 74|) that can amount to pretty large contributions as the order 



of the perturbative expansion increases. Hence, even at high energies, relevant for the Higgs 
decay, the inclusion of these terms to all orders of the perturbative expansion is mandatory, 
albeit a difficult task. It is exactly this issue that singles out the utility of MFAPT because 
within such a perturbative approach, the analytic couplings contain all the aforementioned 
7r 2 terms inherently by construction. 

Our calculation of T(H — > bb) via Rs will be carried out in the MS scheme and will include 
the evolution of both the coupling and the 6-quark mass. Special attention will be given 
below to the origin of the a s corrections in order to distinguish between the loop expansion 
and the loop evolution. In the following, we will compare our results with those obtained in 



the MS scheme in Refs. [63|, [72|, |73j, |74|, where the notation a s = a s /ir (called "couplant") 
was extensively used. For the sake of a better presentation of our results in comparison with 
the existing ones, just mentioned, it is useful to introduce the following abbreviation 

[Ms)ri n = a^\s)^(^Y ^(s). (4.1) 



A. Standard perturbation-theory analysis of Rs 



The Higgs-boson decay into a bottom-antibottom pair can be expressed in QCD by means 
of the correlator 

n(Q 2 ) = (Aixfi ! dxe iqx {Q\ T[ J^(x)J^(0)} |0> 



of two quark scalar (S) currents in terms of the discontinuity of its imaginary part [79 
Rs(s) = Imn(— s — ie)/(27rs), so that the width reads 



T(H bb) 



G f 



Ay/2- 



-M H m 2 b (M K )R s (s = Ml) . 



71 



i.e. 



(4.2) 



Above, Q 2 = —q 2 and J b = ^b^b is the scalar current for bottom quarks with mass nib, 
coupled to the scalar Higgs boson with mass Mh. Direct multi-loop calculations are usu ally 
performed in the Euclidean (spacelike) region for the corresponding Adler function D$ [63j, 
72 . 74| . where QCD perturbation theory works. Hence, we write 



D s (Q 2 ;n 2 ) = 3m 2 b (Q 2 " 



n>l 



(4.3) 



using, as announced above, the notation a s = a s /ir. Connecting to our discussion in Sec. [TT] 
below Eq. (I2.2I) . and adjusting to the scalar case, we now write the results of the 'standard 
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machinery' (see, for instance, 



Rs(s) 



R s (s,s) = 3m b (s) 



ri>\ 



(4.4) 



The coefficients r n contain characteristic L ir 2 terms' due to the (integral) transformation 
R of the powers of the logarithms appearing in D$ (cf. Eq. (I4.3P ). as it was discussed in 
the beginning of Sec. [Ill Recall that in the case of Rs{s), these logarithms stem from two 
different sources: one is the running of a s in Z)g, in analogy to Eq. (12. 2p . the other is related 
to the evolution of ml(Q 2 ). Therefore, the coefficients r n in (I4.4p appear to be related to a) 
the coefficients d n in (12. ip . the latter being directly calculable in the Euclidean space, and 
b) to a combination of the mass anomalous dimension 7$ and the /3-function coefficients bj, 
multiplied by L n 2 powers' 



63j,|72j 
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be substantial, as the following quite recent result, derived in [74 



It turns out that the influence of these ir 2 terms can 

demonstrates: 



[3m 2 ] 1 R s = 1 + 5.667 a. + a 2 [51.57 - 15.63 - N f (1.907 - 0.548 )] 

+ a 3 s [648.7 - 484.6 - N f (63.74 - 37.97 ) + N 2 f (0.929 - 061 



+ a 4 s [9470.8 - 9431.4 - N f 



+N 2 f (54.78 



1454.3 



45.10 ) - N 3 f (0.454 



1233.4 ) 

0.433 )] 



(4.5) 



For emphasis, we have underlined the contributions of those 7r 2 terms which originate from 
the analytic continuation. As one can readily verify, the total amount of these terms is of the 
order of the original coefficients d n , in particular, as regards the coefficient d&. This makes 
it apparent that such terms have to be taken into account in all orders of the perturbative 
expansion. We stress that this is exactly the advantage provided by the analytic machin- 
ery, developed here and in [IJ, and this conceptual advantage arises naturally without any 
additional optimization procedure. Indeed, in FAPT we do not need to expand the renor- 
malization factors into a truncated series of logarithms; instead we can transform them 'as 
a whole' by means of the A M -operation. 

To complete the presentation of the standard analysis, let us display the final result at 
the 0(a% taken from Ref. [jj 



[3ml] 'Ri 



1 + 5.6668 a s + 29.147 a 2 s + 41.758 a 3 s - 825.7 a A s 
1 + 0.2075 + 0.0391 + 0.0020 - 0.00148 . 



(4.6) 
(4.7) 



Note that in Eq. ((¥77]) a s = a 8 (M&) 
boson mass Mh = 120 GeV. 



0.0366 is chosen, which corresponds to the Higgs 



B. FAPT analysis of R s 



We turn now our attention to effects related to the renormalization of the bottom-quark 
mass. For the running mass m^i)(Q 2 ), in the Z-loop approximation, one has the following 
general solution of the RG equation 



m 



(0 



(Q 2 



exp 



a s (Q 2 )/4 



(3(x) 

a s (Q 2 )r f {l) (a s (Q 2 )) 
[a.(A* a )P f(i)M^)) ' 



dx 



(4.8) 
(4.9) 
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where 



Vq = 2— 
bo 



(4.10) 



and the function fm(a s ), given by 



f(l)( c 



exp 



9/4 ( 7m (x) _ Tox 



(4.11) 



accumulates the effects of the second- and higher-loop evolution of m^JQ 2 ) with Q 2 . In the 
one-loop approximation (I = 1), f(i)(a s ) is set by definition equal to unity. On the other 
hand, for / = 2 and I = 3 we obtain 



h ci&o 



f(2)(a s ) = [l + 5i a s ] vi with 8 1 = — 

4o 4 



and 



/(3) (a s ) = [l + 5 1 a s + 5 2 a 2 ] U2 ° exp 



with 



/ 72 _ To 



"21 



v<2x arccos 



-2 6i 



, 7i 7o 

i/1 = 2| A-^ 



1 + Si a s /2 
a/1 + 5ia s + 5 2 > 



72 _ 2 7i , To 



(4.12) 



(4.13a) 



a/4Mo - 6? \ & 2 &i &o 



(4.13b) 



Introducing the RG-invariant quantity rh^, see, e.g., 63. 80j. 

?«(o = ^(o(/^ 2 ) {M/^yp/wM/* 2 ))} 1/2 

one can rewrite Eq. (14.91) in the form 

™ 2 {l) (Q 2 ) = m 2 {l) [a s (Q 2 )Y° h)MQ 2 )) , 
where the expansion of f(n(x) at the three-loop order is given by 



(4.14) 



(4.15) 



/(0( ( 



1 _l fo i /^7i To , 



2 °1 



s 16b 2 



To Ti \ , M2 / T2 To 



b h 



b\ 



+ o (o2) 



&2 &0 

(4.16) 



which is in one-to-one correspondence with Eq. (15) found by Chetyrkin in [811] • [Note, 
however, that Chetyrkin expands instead the expression \J f(i)(x) = c(x) and uses different 
normalizations; viz., j3 n = b n j (4™6 ) an d Tn = Tn/ (4"&o)-] Keeping in mind the main purpose 
of our task, namely, the sequential "analytization" of D, we rewrite the RG Eq. (14.151) in 
the form of a power series to get 



m 2 (l) (Q 2 ) 



m 



2\\ u 



(a s (Q 2 )) 



i + E e 2 «Q 2 )Y 



m>l 



(4.17) 
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where the coefficients em depend implicitly on the RG parameters via Eq. (14.111) . For the 
simplest case of a two-loop running, they can be written down explicitly: 



r(m + i)r(i/i -m + iy h 



(4.18) 



Substituting Eq. (fl~TTl) for m 2 {l) {Q 2 ) into D S (Q 2 ) = D S (Q 2 ] Q 2 ) in Eq. (TOD, one finds 

l oo 



n>l 



m>l 



with 



min[Z,m— 1] 



A« = e (l) + V d t e 



(0 

m—k 



(4.20) 



fc>i 



The effects of mass evolution of the higher orders are collected in the third term on the RHS 
of Eq. (I4.19p and have been purportedly separated from the original series expansion of D 
(truncated at n — I), the latter being represented by the second term on the RHS of Eq. 
fl4.19p . In practice, for Q > 2 GeV, i.e., for ot s < 0.4, the truncation at m = I + 4 of the 
summation (14.171) produces a truncation error smaller than 0.01%. 

To obtain R^ FAPT , we recall the action of the Am operation in FAPT, as described in 
Sec. HIl and consider the map of the quantity Ds\Q 2 ) in Eq. (14.191) onto the Minkowski 
region. Following the "analytization" procedure illustrated in Fig. [2j we then obtain 



R 



(Z)MFAPT 



Am[4°] 



3 rh 



(I) 



n>l 



m>l 



(4.2i; 



where the superscript I denotes the loop order of the evolution and fixes at the same time 
the order of the perturbative expansion of the Dg-function. The above expression contains, 
by means of the coefficients A„ ( ej^ ) and the couplings <x^ +vo , all RG terms contributing 
to this order, while the resummed n 2 terms are integral parts of the analytic couplings by 
construction. 



C. Comparison of different perturbative approaches to obtain Rs 

We list below the results obtained for the quantity R$ using different methods. 

• Broadhurst, Kataev, and Maxwell (BKM) [63j utilized within the so-called "naive 
non-Abelianization" (NNA) approach an optimized power-series expansion, based on 
the "contour integration" technique, to compute Rs- Their estimate, with one- loop 



19 



running of a s and setting ai ^ = a s , reads (see Section 3.3 in [63j) 



]\3 



R 



(i=l)BKM 



.4 



BKM / 



3m 2 (l=1) (a s f° 



1 



BKM, 



b 7i5 n 
n + UQ — 1. 



b 7r a t 



n>l 

Sn/2 



,n-i ■ I K , I b na s 
a, sin o n arctan 



(4.22) 

(4.23) 
(4.24) 



These new couplings A^ KM (a s ) and the whole result are closely related to our analytic 
approach at the one-loop level, as we will show shortly. 

Baikov, Chetyrkin, and Kiihn (BChK) [74j have derived within the standard pertur- 
bative QCD at the 0(a^), c.f. Eq. (14.51) ). the following expression: 



£(*=4)BChK = 3 ^2 (s) G=4) 



n>l 



(4.25) 



Consider now the MFAPT equation (14.211) and recast it in the form 

l 1+4 



R 



(i)MFAPT 



3 rh 



(i) 



m>l 



3 rh 



(i) 



J,r rf fl (0 +VA w a (0 

n>\ 



"0 



1 (0 

n+i't) 



4 \ 4 

Ve ( V° +Ve (0 a 

/ j °m u n+m+uo J ~ / j °m " 



(0 

m+VQ 



m>l 



(4.26) 
(4.27) 



where we have truncated the "evolution series" at m — I + 4, adopting the empirical 
recipe discussed after Eq. ( I4.20p . Pay attention that the terms a^ +i , contain — by 
means of the index i> Q — all 70 and 60 terms, and also all tc 2 terms, while the con- 
tributions of the higher-loop RG-dependent parts are accumulated in the coefficients 
AS— see Eqs. (Q5j) . (14161) and also (l4~T9j) . (ET20]) — in terms of the RG-parameters 7; 
and by 



Equation (I4.26P can be considered as a generalization of the BKM "contour-improved" ex- 
pansion in Eqs. (I4.22j) (I4.24|) . To see this, one should reduce the above expression to the 
one-loop case, by setting Am — > and a$ 



, with the aim to reproduce the "contour- 



improved" effective coupling in Eq. (I4.22p . In fact, taking into account the explicit expression 
(13. 2p for ai , one can readily find the relation 



bo 



n+uo 
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(i) 

n+uo 



(a s y°A n (a s ) 



(4.28) 



3 The couplings a s and a„ are understood to be functions of s, i.e., a s (s) and a^s). Note in this context 
that in the original paper of Ref . [63[ , the authors used the coefficients d^ NA , which are the "all-order" 
coefficients of the expansion of the Adler function in the Euclidean region, estimated, however, through 
the NNA procedure. 
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recalling (I4.23P and (14.241) . This relation establishes the equivalence between the "contour- 
improved" effective coupling and the timelike analytic coupling in the one-loop approx- 
imation of MFAPT. The only difference between ^- 1 ) BKM anc [ ^(J— i)mfaft s ^ emg f rom 
summation in Eq. (I4.22p . which extends to all known ^-coefficients, i.e., also to those up to 
n = 4. 

In Fig. [TJ we display the final results for the quantity ^W MFAPT j n Eq.( |4.26p vs. the Higgs 
mass M H evaluating it for the cases of two (/ = 2) and three (/ = 3) loops, with the goal 
to illustrate the effects of the "analytization" procedure in comparison with the standard 
approach. The solid blue line in both panels of this figure shows the prediction obtained 
with MFAPT and a fixed number of active flavors Nf = 5. One appreciates that this curve 
lies only slightly (about 2% for I = 3) above the standard result i?g _4 ' )BChK , illustrated by 
the dashed red line. This is mainly due to the somewhat larger values of the perturbative 
coefficients d n within MFAPT relative to the standard ones, r n . At the same time, the dotted 
green line, which corresponds to i?g 1 , Eq. (4.22), turns out to lie lower than the BChK 
prediction by about 7 to 8%. Note that if one sets Ayy,=5 = 221 MeV, as dictated by the one- 
loop perturbative QCD analysis of the DIS data [82], then the BKM curve, corresponding 
to this value, will end up to lie higher than the BHhK prediction by 8%. 



34 



32 



30 



28 - 



26 - 



- R s (Mn) [GeV 2 ] 






M H [Gevf~* 



34 



32 



30 



28 



26 




60 



80 



100 



120 



140 



60 



80 



M H [GeV] 

100 120 



140 



FIG. 7: Illustration of the calculation of the perturbative series of the quantity Rs(M^) in different 
approaches within the MS scheme: standard perturbative QCD 73|, |74| at the loop level I = 4 
(dashed red line, Atv /= 5 = 231 MeV); BKM estimates, by taking into account the 0((a s ) u ° A4(a s ))- 
terms, (dotted green line, Aat /=5 = 111 MeV); and finally MFAPT from (TOTl) for Nf = 5 
(solid blue line), displayed for Z = 2 (left panel, Aat /= 5 = 263 MeV) and / = 3 (right panel, 
Ajv /= 5 = 261 MeV). All Ajv /= 5 scales are defined in such a way as to reproduce the corresponding 
couplings with the value 0.120 at s 



m 



In TableUJwe show for the quantity i?s( s ) a comparison of the perturbative series evaluated 
in different expansion approaches, discussed in the literature and in this work. One may 
conclude that the standard perturbation-theory series and MFAPT show a similar behavior, 
starting with the two-loop running. Because the coefficients d n of MFAPT are larger than 
the standard ones, r n , due to the subtraction of the it 2 terms in the latter, the fixed-order 
MFAPT result for the quantity Rs{s) appears to be slightly larger than that of the standard 
perturbative QCD approach. However, the coupling parameters of MFAPT, a u , contain the 
resummed contribution of an infinite series of 7r 2 -terms that renders them ultimately smaller 
(for the same value of Aqcd) than the corresponding powers of the standard coupling, as 
it can be seen from Table [Til In order to demonstrate the combined effect of evolving the 
coefficients and the heavy quark mass at different orders of the perturbative expansion, we 
define the following factors for, respectively, the standard perturbation theory (PT) and for 



21 



TABLE I: Comparison of the perturbative-series convergence in terms of the associated pertur- 
bative orders, using for the Higgs boson mass the value Mh = 120 GeV and employing different 
expansions within the MS scheme: (i) Standard perturbative QCD [73|, [7J], (ii) BKM 63], and 



(iii) MFAPT at the two- and the three-loop level of evolution with the number of flavors fixed to 
N f = 5 and s = [120 GeV] 2 . 



Expansion approach used 


Rs(s) 


0(1) 


0(a s ) 


0(al) 




o{<4) 


standard QCD [74] 


27.44 GeV 2 


80.2% 


16.6% 


3.1% 


0.2% 


-0.1% 


MFAPT 


Ms) 




o(«SL) 




o(4L) 




BKM [63] a (MFAPT for I = 1) 


31.89 GeV 2 


74.5% 


17.7% 


5.3% 


1.8% 


0.7% 


MFAPT for / = 2, this work 


27.59 GeV 2 


79.5% 


16.2% 


4.3% 






MFAPT for I = 3, this work 


28.07 GeV 2 


78.5% 


16.1% 


4.2% 


1.2% 





"Note that in this order of expansion, the analytic couplings in MFAPT, a v , and those in the "contour- 
improved" approach coincide, though the associated coefficients d n and d^ NA , respectively, are different. 

MFAPT: 



(4.29) 
(4.30) 



m>l 



Then Eqs. (14.251) and (14.271) can be rewritten as 



n(')PT _ q - 2 



f MFAPT = 3m 2 



(i) 



n>l 
I 



n>l 



(4.31) 
(4.32) 



In order to understand the effect of "analytization" , we compare the values of £n+ v As) with 



those of As), and display their ratio in Table HT1 Using the entries of this Table, one 



TABLE II: Comparison of the combined evolution effects on the ratio of MFAPT to the standard 
perturbation theory for the three-loop evolution case and fixing the number of flavors to Nf = 5. 
For the sake of completeness, we also show the ratios of the associated perturbative coefficients d n 
and r n and the approximate relative weights uj n = r n E^ UQ / fiE^ of the corresponding 



i=0,...,3 

contributions. 





n = 


n = 1 


n = 2 


n = 3 


r<?(3) /p(3) i 

[°n+u 1 ^n+u \ s =(i20 GeV) 2 
dn/ fn 

relative weight oj n 


1.00 
1.00 
0.801 


0.98 
1.00 
0.166 


0.95 
1.44 
0.031 


0.91 

8.46 
0.002 
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can easily estimate the relative enhancement of ^ 3 ) MFAPT w j^ n reS p ec t to R^ PT : 

5(3)MFAPT 3 , ~(3) / x 

s ~„ ^ = V u n — n + uo{ f = 0.801 + 0.163 + 0.042 + 0.015 = 1.021 , (4.33) 
where u n = r n E®J ^ r t E^ Uo . (4.34) 

i=0,...,3 

From this equation in conjunction with the values of u n (presented in Table HT]l . we see 
that the largest enhancement is provided for n = 3 due to d n /r n = 8.46, amounting to a 1.5% 
contribution to ^ 3 ) MFAPT ou t Q f 2.1% in total. Evolution, which potentially could reduce 
this enhancement, because of the inclusion into the analytic coupling of the resummed tt 2 
terms, notably, £n+ Uo /En+v — 0.91, is too small to counterbalance it. 



V. SUMMARY AND CONCLUSIONS 

This report has focused on the implementation of analyticity of QCD amplitudes both in 
the Euclidean and in the Minkowski space, the theoretical basis being provided by dispersion 
relations (to ensure causality) in conjunction with the renormalization group. The goal of 
the work was to elevate Analytic Perturbation Theory — initiated by Shirkov and Solovtsov 
-to a calculational paradigm for perturbative QCD applications, capable of providing 



singularity-free expressions for any real power of the strong coupling. Following the rationale 
that all quantities that ma y c ontribute to the spectral density should be included into the 



analytization" procedure [6ll . [62] , we have expanded the previously developed Fractional 



Analytic Perturbation Theory [1] from the spacelike to the timelike region. 

The core issues of the investigation can be summarized as follows, starting with the 
benchmarks of the analytic couplings Av\l) (spacelike) and 21? (L s ) (timelike) of (M)FAPT 
for an arbitrary real index v: 

• j*4?(L), 2r?(L s ) j is analytic in L (respectively, L s ) and also in the index i/GZ; for 
I = 1, see Eqs. (13.91) and (13 . lOj) : for / = 2, see Sec. Ill C and also [27| (for integer 
indices). This implies the existence of any index derivative — — |.4.?(L), 21? (L s )}. 

• A^(-oo) = 2l#(-oo) = 5 hm form6N0 and A2\-oq) = S^°(— oo) = for v > 1. 
This implies the existence of a universal IR fixed point. 

• 21? (L s — > oo) — > Ai (L — > oo) — > clq\(L — > oo) have the correct UV asymptotics, 
dictated by the asymptotic freedom of QCD. 

• The effect of the "distorted- mirror symmetry" with respect to L and L s , is exhibited 
analytically in terms of Av\l) and 21? (L s ) for any real v. 

• The "analytization" of expressions, like 4l(L) = A E [(a (l) (L)) u L m ], £$ m (L,) = 
Am [{ a (i)(L)) u £ m ] , which amount to 'evolution logarithmic factors' in /-order QCD 
perturbation theory, typical examples being logarithms of the factorization scale (ap- 
pearing beyond the leading-order expansion), and significantly affecting the precision of 
perturbative calculations, becomes possible. Moreover, the insensitivity to the choice 
of the factorization scale [67[ and a diminished dependence on the adopted renormal- 
ization scheme and scale setting has been achieved (36| . 
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The following topics have been given particular attention: 



1. We have presented a brief historical review of the "analytization" technology from early 
attempts in QED up to the most recent developments in QCD perturbation theory and 
phenomenology (Introduction), highlighting the strategy for resolving the conflict with 
analyticity in Q 2 of QCD "observables" in the perturbative regime. 

2. In sections HT1 and HTT1 we have worked out in detail the theoretical framework to deal 
with analytic versions of the strong coupling and its powers for any fractional (real) 
power, encompassing the spacelike and also the timelike region. 

3. We have given closed- form expressions for the analytic-coupling images at the one- loop 
level both in the Euclidean, as well as in the Minkowski space. We have also provided 
explicit but approximate expressions at the two-loop level, which show excellent agree- 
ment with the exact but numerical results in terms of the Lambert function, found 



before by Magradze 25] — see Fig. Three-loop (approximate) results have also been 



included and the fast convergence of the non-power-series expansion has been proved. 

The major characteristics (advantages) of this machinery, studied in Sec. 1111} are: 
(i) a reduced uncertainty level of only a few per mil in a wide range of momentum 
(or energy) values, starting just above 1 GeV, (ii) the high quality of the non-power- 
series expansion, evidenced by providing trustworthy analytic expressions for the strong 
coupling and its powers at the one, two- and three-loop level that are singularity free 
in the spacelike region and resum the 7r 2 terms (induced by analytic continuation) to 
all orders in the timelike domain. 

The relevance of the developed framework for practical purposes in the Minkowski 
space has been effected in Sec. [TV] by applying it to the decay of a scalar Higgs boson 
to a bb pair, this task serving as a Proof-of-the-Concept-calculation. Specifically, we 
estimated the quantity i?s at the four- loop level of the perturbative expansion (i.e., 
up to the coefficient d 3 ), having recourse to the coefficients and anomalous dimensions 



calculated by Chetyrkin and collaborators in [72], [73J . The main advantage of MFAPT 
here is that the coupling parameters 2l„ inherently include the resummed contribution 
of an infinite series of those 7r 2 -terms originating from the analytic continuation from 
the Euclidean to the Minkowski space. Technically, we could have also included the 
next correction, associated with the coefficient d^, but this would be not worth the 
effort, given that the expected contribution would be about 0.5%. 

In conclusion, we think that the analytic approach to QCD perturbation theory in the 
form advocated here has indeed been rather successful both in the spacelike and in the 
timelike region for the particular processes we have discussed in this work and elsewhere [U, 
36, 63]. The developed full-fledged machinery may improve our understanding of key issues 



of QCD reactions from the poorly understood low-momentum (spacelike) domain, where the 
Landau singularities are a serious obstacle in standard perturbative QCD (based on power- 
series expansions) up to the high-energy (timelike) region of several GeV. In conjunction with 
high-loop calculations in the latter case, the 7r 2 terms, entailed by analytic continuation, 
can significantly change the result in higher orders. A greater challenge, however, is to 
include into the approach gluon resummation and power corrections analytically, beyond 
the attempts in [6l| and [lj, with the aim to unravel the analytic structure of the involved 
amplitudes in the whole momentum (energy) range. Work in this direction is currently in 
progress. 
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APPENDIX A: TWO-LOOP RENORMALIZATION-GROUP SOLUTIONS FOR 
THE COUPLING 



1. The expansion of the /3-function is given by the RHS of the equation 

£(£) ->(£)- (£)'-•• (£)*-■■■ 

where L = ln(/i 2 /A 2 ) and 



(Al) 



b 
b 2 



— C A - - T R N f 
3 3 



6i = jC%- (4C F + jC A )T R N f - 



2857 
~54~ 



C% + 2C*T R N f 

2 



205 



C F C A T R N 



f 



1415 



44 



ClT K N f + -C F (T R N 



+ 1 -§C A (T R N f) , 



(A2) 



with C F = (A 2 - 1) /2A C = 4/3, C A = N c = 3, T R = 1/2, and N f denoting the number of 
active flavors. The corresponding three- loop RG equation for the coupling a = b a /{An) is 



da 



(3) 



dL 



-0(3) [1 + ci op) + c 2 af 3) ] with c x = ^, c 2 = 



(A3) 



The solution of this RG equation at the two-loop level (c 2 = 0) assumes the form 



1 

1(2) 



+ c\ In 



1(2) 



1 + Ci<Z(2) 



L. 



(A4) 



The exact solution of Eq. flA4j) can be expressed in terms of the Lambert function W{z) 
[83| (see also 0, 23]) defined by 



z = W{z) e w{z) 



(A5) 



This solution has the form 



a (2 )(L) 



Cl l + W_i(«(L)) 



(A6) 
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where z(L) = (1/ci) exp (—1 + in — L/ci) and the branches of the multivalued function W 
are denoted by Wk, k = 0, ±1, .... A review of the properties of this special function can be 
found in j83|; see also [22| for a special emphasis on the problem considered here. 

2. The expansion of the solution of Eq. (IA30 . i.e., a( 3 )(L), in terms of a = l/L, while 
retaining terms of order a 4 , yields 



1(3) 



+ a 2 c\ In a(i) + a 3 [c\ (in 2 a + In a — l) + c 2 ] 



(In a + — In a — 2 In a — — J + 3c2Ci In a 



+ a 

+ (a 5 In 4 a) . 

For any power v of the coupling we have the final three-loop expression 

y + l 



(A7) 



(<J(3)) 



a + Ci^a ma + 



2! 



+ c 3 z/a"+ 3 



(z/+ l)(i/ + 2) 



In a + In a — 1 ) + C2 

2^ + 3 



z/ a 



1^+2 



In 3 a H — In 2 a - (v + 1) In a - 



+ Cic 2 v(v + 2) a" +3 lna 



+ c 4 z/a^+ 4 



(i/ + + 2)(i/ + 3) 4 3 z/ 2 + 12 i/ + 11 



4! 



In 4 a + 



z/ 2 + 2z/ 2 3z/ + 5, 

In a In a + 



6 

3z/ + 4 



In 3 a 



+ c 2 v a 



2 

+ 4)c 2 



3^ + 7 2 



+ c 2 {y + 2) 



6 

z/ + 3 



In a + In a — 1 



+ (a" +5 In 5 a) . 



(A* 



Note that the two-loop expansion can be immediately obtained from Eq. (IA8I) by setting 
c 2 = 0. 

3. For completeness, we also present here the expansion of the product [a (2)] £ i n terms 
of a that can be obtained (in the two-loop approximation) from Eq. (1A4j) : 



(a( 2 ))"L= (a {2) ) u 1 + (a (2 ))"ci In ^ — ) • 
Expanding the logarithmic term In (l + Cid( 2 )) , while retaining terms of order 



(A9) 



(2) ' 



a p {2) ln(a {2) ), ajg , ajj , we get 



(a (2) )^L = (a(2)) !/ X + Ci (a (2 ))^ In a (2 ) - c 2 a^ 4 + ya^ 2 + O (a^ '' 



'(2) 



(A10) 



and, finally, expanding the coupling ar 2 ) in terms of a = am, 



c\ v a u In a + c 2 z/ a 
+ (a^ +3 In 3 a) 



H-i 



^+1,2 



ln 2 (a) + ln(a) - 1 



(AH) 
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in accordance with Eq. (1A8I) . An analogous expression can be constructed for the three-loop 

case: 

V + l 



Ci v a u In a + 



2! 



In a + hi a — 1 + c 2 



i/ a 



f+1 



[v + l)(v + 2) ]3 2i/ + 3 . 2 
— In a H hr a 



Cic 2 z/(z/ + 2) a 1 ^ In a 



3! 

u+2 ■ 



Ina-- 



{v + l)(u + 2){u + 3) , 4 3 z/ 2 + 12 i/ + 11 , 3 

71 In a H In a 

4! 6 

v 2 + 2v , 3i/ + 5, 3z/ + 4 
In a — In a + 



,W-3 



+ c 2 v a 
+ O (a^ 4 In 5 a) . 



z/ + 4)c 2 ^±^ + c 2 (z/ + 2) 



6 



6 

u + 3 



In a + In a — 1 



(A12) 



4. We consider here the solution of the renormalization-group equation in the so-called 
Pade-modification of the three- loop approximation in QCD |24j, |26|, where the /3-function of 
Eq. (IAip is given by 



6i (a a /(47r)) 



6 (1 - 6 2 a s /(46i7r)) 
The corresponding three-loop RG equation (1A3|) is modified to read 



da 



(3-P) 



l (3-P) 



1 + 



Cl Q(3-P) 



1 - (C 2 /Cl) 0(3-P). 



The solution of this RG equation assumes the form 



1 



«(3-P) 



Ci In 



0(3-P) 



(A13) 



(A14) 



(A15) 



1 + (ci -c 2 /ci)a(3_ P) _ 

which is very similar to Eq. flA4j) . For this reason, the exact solution of Eq. (1A15j) can also 
be given in terms of the Lambert function W(z). The solution is 



Ci 1 — 



c 2 /cl + W^{z^){L)) 



(A16) 



where z^ 3 ~ F \L) = (l/cj.) exp [—1 + in + c 2 /c\ — L/ci). The relative accuracy of this solution, 
as compared with the original Eq. flA3j) solved numerically, is better than 1% for L > 7 (and 
better than 0.5% for L > 9). 

APPENDIX B: SPECTRAL DENSITY AT HIGHER-LOOP LEVEL 



1. We consider here the spectral density p u (a) beyond the one-loop approximation. At 
the /-loop level, Pu\cr) can always be presented in the same form as for the leading order, 
given in Eq. (12. 12ft . i.e., 

_ 1 sin[z/ <fi(i)(a)) 



7T W 



Bl 
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FIG. 8: The dashed red line corresponds on the left panel to the first running-coupling iteration, entering 
the spectral density, p^' 1 *' 1 (L a ), see Eqs. (B12)-(B14) in and, on the right panel, to the second running- 
coupling iteration, p„l~l t 2 {L a ), while the solid blue line on both panels corresponds to the exact spectral 

For the sake of a better comparison, the displayed region of L a is chosen in 



density p\ (L a ), see Eq. 
the vicinity of the maximum difference of these two curves. 



where now the phase (pn\ and the radial part Rn\ have a multi-loop content. At the two-loop 
level, one has to consider the imaginary part of the Lambert function W-\ (cf. Eq. (1A6j) ); 
see later discussion. We consider first the known second-order iterative solution of Eq. (1A4I) 
that provides us, with sufficient accuracy, the following result: 



L + ci In [L + ci + ci In (L + a)] 



For the approximate solution a^ 2 f, we have 



[L a + Cl lni?(6r)] 2 + 
L a + c\ lni?(<j 



?T + c\ arccos 



L a + ci + ci In r(a) 
W) 



arccos 



i?f 2 f(a) 



where 



R(a) 
r(a) 



\L<j + ci + ci lnr(a)] 2 + [it + c l( f>(a)]' 



[L a + Cl f + 7T 2 



V) 



arccos 



r(a) 



with L a = ln(cr/A 2 ). The spectral density 



a=2)it- 2 , , i 8111 



71 



Dit-2, 
-"■(2) 1 



(J 



(B2) 



(B3) 
(B4) 



(B5) 
(B6) 



(B7) 



with the phase ip^(a) and the radial part R^(a) from Eqs. (1B3|) - (1B6I) . appears to be very 

(2) 

close to the exact, but numerical result for p\ (cr), based on the Lambert function W-\ — see, 



e.g. 
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Specifically, one can use the symbolic program Mathematical or Maple 7, which both 
recognize the Lambert function, to carry out these integrations using the spectral density 
|2T 



PmV) = -Im 



7T 



-1 



Ld (i + w_i(^)) 



(B8) 



where z a = c x 1 exp (— L a jc\ — 1 We present the approximate expression p^ 11 ' 2 (L a ) 

(dashed red line) in comparison with the exact expression pf\L a ) (solid blue line) in Fig. El 
2. It is worth noting here that explicit expressions for the analytic images of the integer 
powers of the couplings in the Minkowski region 



< 2) (Ls 



p®{a)dL a 



(B9) 



have been obtained before by Magradze j27J, using the properties of the Lambert function; 
notably, 



2lS 2) (L 

f(2) 



1 - - Im pnWiCz,)] 

7T 



1 



%%>[L a } = Im 

CxTT 



111 



W 1 (z s 



1 + W x (z s 



1 

Ci 



for n > 1 . 



(BIO) 
(Bll) 
(B12) 



3. For the Pade- modification of the three-loop approximation the corresponding spectral 
densities are defined through Eq. (IA160 



7T 



•1 



cx (l-c 2 /c\ + W„x{z^)(L a ))) 



(B13) 



The following explicit expression for the analytic image of the coupling in the Minkowski 
region 



n { cf-c 2 



lnWx(z s 



° 2 T 
+ -r, Im 

c{ - c 2 



In f l-% + W x (z a 



(B14) 



with z s = z( 3 ~ p \L s ) has been obtained before by Magradze [23(. The relative accuracy of 
this solution, as compared with the numerical integration of the original, non-Pade-modified, 
spectral density p[ (cr) is better than 0.25% for L > 2. 

4. A useful formula to relate arbitrary powers of logarithms by means of the dispersion 
relation (12.61) has been presented in (63|; viz., 



Q 2 



ds 



(s + Q 2 ) 2 \s 



7l5 ( jl 



sin(7r5) \Q 2 



(B15) 



In versions 3, 4, and 5 of Mathematica the function Wk(z) is denoted by the name ProductLog[fc, z]. 
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APPENDIX C: "ANALYTIZATION" OF POWERS OF THE COUPLING MUL- 
TIPLIED BY LOGARITHMS 



1. Here we derive the "analytization" of powers, (a(2)Y ', and more complicated expres- 
sions that contain powers of the running coupling multiplied by logarithms, making use of 
the property 



[a I/ ln(a)] E 



d 

dv 



def d 

dv 



(CI) 



supplemented by Eqs. (IA7[) and (|A8j) . In this way we obtain (ci = c 2 = b 2 /bl) 
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2. The images of the coupling accompanied by logarithms of the momentum, i.e 
(ci(2)) A hi accordance with Eq. (lAlOj) . are 

+ O (A?l 3 ) . (C3) 
Following Eq. ( 1A12[) . this leads to 
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APPENDIX D: REPRESENTATION OF THE D-FUNCTION 



1. The first three coefficients d x ,d 2 , d% for the D$ function of two scalar quark currents, 

)2\ \ n 



D S (Q' 2 ) = 3m 2 b (Q 2 " 



n>0 V 7 



(Dl) 



30 



have been calculated in [721 and read 
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The result for c?4 was obtained quite recently in [74 
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2. The coefficients ji determine the expansion of the quark- mass anomalous dimension 
in analogous manner as Eq. flAlj) does with respect to the expansion of the /3-function; viz., 
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with the following expressions, see [81 
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where £(z/) denotes the Riemann £ function. 
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